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Abstract. The standard central limit theorem plays a fundamental role in 
Boltzmann-Gibbs statistical mechanics. This important physical theory has 
been generalized |T in 1988 by using the entropy S q = t\ p - (with qeTZ) 
instead of its particular BG case Si = Sbg = — Y2% Pi m Pi- The theory which 
emerges is usually referred to as nonextensive statistical mechanics and re- 
covers the standard theory for q = 1. During the last two decades, this q- 
generalized statistical mechanics has been successfully applied to a consid- 
erable amount of physically interesting complex phenomena. A conjecture[2] 
and numerical indications available in the literature have been, for a few years, 
suggesting the possibility of (/-versions of the standard central limit theorem 
by allowing the random variables that are being summed to be strongly cor- 
related in some special manner, the case q — 1 corresponding to standard 
probabilistic independence. This is what we prove in the present paper for 
1 < q < 3. The attractor, in the usual sense of a central limit theorem, is 
given by a distribution of the form p(x) — C q [l — (1 — q)f3x 2 ] 1 ^ 1 ~ q ^ with 
P > 0, and normalizing constant C q . These distributions, sometimes referred 
to as g-Gaussians, are known to make, under appropriate constraints, ex- 
tremal the functional S q (in its continuous version). Their q = 1 and q = 2 
particular cases recover respectively Gaussian and Cauchy distributions. 
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1. INTRODUCTION 

Limit theorems and, particularly, the central limit theorems (CLT), surely are 
among the most important theorems in probability theory and statistics. They 
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play an essential role in various applied sciences, including statistical mechanics. 
Historically A. de Moivre, P.S. de Laplace, S.D. Poisson and C.F. Gauss have first 
shown that Gaussian is the attractor of independent systems with a finite second 
variance. Chebyshev, Markov, Liapounov, Feller, Lindeberg, Levy have contributed 
essentially to the development of the central limit theorem. 

It is well known in the classical Boltzmann-Gibbs (BG) statistical mechanics 
that the Gaussian maximizes, under appropriate constraints, the Boltzmann-Gibbs 
entropy Sbg = — ^iPi^Pi (Sbg = — J dx p(x) hi p(x) in its continuous form). 
The g-generalization of the classic entropy introduced in pQ as the basis for gener- 



alizing the BG theory, and given by S q — 1 ^ i j Pi (S q 



1-Jdx[p(x)}1 /(q-1) 



1- 

in its continuous form), with q £ 1Z and Si = Sbg, reaches its maximum at the 
distributions usually referred to as g-Gaussians (see [3]). This fact, and a number 
of conjectures [2], numerical indications [3], and the content of some other studies 
[U El [5 . 6 suggest the existence of a g-analog of the CLT as well. In this paper we 
prove a generalization of the CLT for 1 < q < 3. The case q < 1 requires essentially 
different technique, therefore we leave it for a separate paper. 

In the classical CLT, the random variables are required to be independent. 
Central limit theorems were established for weakly dependent random variables 
also. An introduction to this area can be found in [TJ [HI El El [HI H2] (see also 
references therein) , where different types of dependence are considered, as well as 
the history of the developments. The CLT does not hold if correlation between 
far-ranging random variables is not neglectable (see [13]). Nonextensive statistical 
mechanics deals with strongly correlated random variables, whose correlation does 
not rapidly decrease with increasing 'distance' between random variables localized 
or moving in some geometrical lattice (or continuous space) on which a 'distance' 
can be defined. This type of correlation is sometimes referred to as global correlation 
(see [IT] for more details). 

In general, g-CLT is untractablc if we rely on the classic algebra. However, 
nonextensive statistical mechanics uses constructions based on a special algebra 
sometimes referred to as g-algebra, which depends on parameter q. We show that, 
in the framework of g-algebra, the corresponding ^-generalization of the central 
limit theorem becomes possible and relatively simple. 

The theorems obtained in this paper are represented as a series of theorems 
depending on the type of correlations. We will consider three types of correlation. 
An important distinction from the classic CLT is the fact that a complete formula- 
tion of g-CLT is not possible within only one given q. The parameter q is connected 
with two other numbers, q* = z~ x (q) andg* = z(q), where z(s) — (l+s)/(3— s). We 
will see that identifies an attractor, while q* yields the scaling rate. In general, 
the present ^-generalization of the CLT is connected with a triplet (qk-i, qk, Qk+i) 
determined by a given q G [1,2). For systems having correlation identified by qk, 
the index qk-i determines the attracting q-Gaussian, while the index qk+i indi- 
cates the scaling rate. Note that, if q = 1, then the entire family of theorems 
reduces to one element, thus recovering the classic CLT. 
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The paper is organized as follows. In Section 2 we recall the basics of q- 
algebra, definitions of q-exponential and q-logarithm. Then we introduce a trans- 
form F q and study its basic properties. For q ^ 1 F q is not a linear operator. Note 
that F q is linear only if q = 1 and in this case coincides with the classic Fourier 
transform. Lemma 12.51 implies that F q is invertible in the class of g-Gaussians. 
An important property of F q is that it maps a g-Gaussian to a g*-Gaussian (with 
a constant factor). In Section 3 we prove the main result of this paper, i.e., the 
g-version of the CLT. We introduce the notion of ^-independent random variables 
(three types), which classify correlated random variables. Only in the case q = 1 
the correlation disappears, thus recovering the classic notion of independence of 
random variables. 



2. g- ALGEBRA and g-FOURIER TRANSFORM 
2.1. g-sum and g-product 

The basic operations of the g-algebra appear naturally in nonextensive statistical 
mechanics. It is well known that, if A and B are two independent subsystems, then 
the total BG entropy satisfies the additivity property 

S BG (A + B) = S BG (A) + Sbg(B). 

Additivity is not preserved for q ^ 1. Indeed, we easily verify [UGH [2] 

S q {A + B) = S q {A) + S g (B) + (1 - g) S q (A) S q {B). 

Introduce the g-sum of two numbers x and y by the formula 

x® q y = x + y+(l- q)xy. 

Then, obviously, S q (A + B) = S q (A) ® q S q (B). It is readily seen that the g-sum 
is commutative, associative, recovers the usual summing operation if g = 1 (i.e. 
x 0i y = x + y), and preserves as the neutral element (i.e. x (B q = x). By 

inversion, we can define the q-subtraction as x Q q y = ■ Further, the q- 

i 

product for x, y is defined by the binary relation x ® q y — [x 1 ^ 9 + y x ~ q — 1] jj_ ^ . 
The symbol [x) + means that [x] + = x, if x > 0, and [x] + = 0, if x < 0. Also this 
operation is commutative, associative, recovers the usual product when g = 1 (i.e. 

x (£>i y = xy), preserves 1 as the unity (i.e. x <giq 1 = x). Again, by inversion, a 

i 

q-division can be defined: x Q) q y = [x 1 " 9 — y x ~ q + l]^" ■ Note that, for g ^ 1, 
division by zero is allowed. As we will see below, the g-sum and the g-product 
are connected each other through the g-exponential, generalizing the fundamental 
property of the exponential function, e x+v = e x e y (or hi(xy) — In a; + lny in terms 
of logarithm) . 
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2.2. ^-exponential and q-logarithm 

The g-analysis relies essentially on the analogs of exponential and logarithmic 
functions, which are called ^-exponential and q- logarithm [18] . In the mathematical 
literature there are other generalizations of the classic exponential distinct from 
the g-exponential used in the present paper. These generalizations were introduced 
by Euler [T2], Jackson [T5], and others. See [T5] for details. 

The q- exponential and q-logarithm, which are denoted by e x q and ln q x, are 



and hi q x — 



(x > 0). 



= elel and e x+y = 



'9 9 



respectively defined as e x = [1 + (1 — q)x] + " 

For the q-exponential, the relations e\ 
hold true. These relations can be rewritten equivalently as follows: ln ? (x ® q y) 
hi q x + hi q y, and ln q (xy) = ln q x(S q ln q jQ. The (/-exponential and (/-logarithm have 
asymptotics e x = l+x+^x 2 +o(x 2 ), x — * and ln g (l+a;) = x— ^x 2 +o(x 2 ) 7 x — > 0. 
For q ^ 1, we can define e q v , where i is the imaginary unit and y is real, as the 
principal value of [1 + i(l — q)y] T -* , namely 

,-, .-. 1 i arc tan \ (q— l)v] 

e\y = [1 + (1 - q) 2 y 2 ] — e ^ , q + 1. 

Analogously, for z = x + iy we can define e q — e x+zy , which is equal to e x ® q e l q v in 
accordance with the property mentioned above. Note that, if q < 1, then for real y, 
\ef\ > 1 and \e q v \ - (1 + y 2 )^) , x -> oo. Similarly, if q > 1, then < \ef \ < 1 
and \ef\ -> if |y| -> oo. 



2.3. (/-Gaussian 

Let j3 be a positive number. We call the function 



VP - 



(ix 1 



(2.1) 



a q- Gaussian; C q is the normalizing constant, namely C g = j l _ 
to verify that 



e g x dx.lt is easy 



== [„ ' (COS t) 1-9 dt = r 



Cq = { 



7=t/ (i + r)'- 1 ^y = 



g-ir 



— oo < q < 1, 

9 = 1, 
1< g < 3. 



(2-2) 



For g < 1, the support of G q ((3\x) is compact since this density vanishes for 
1 03 j > — q)/3. Notice also that, for q < 5/3 (5/3 < q < 3), the variance is 

finite (diverges). Finally, we can easily check that there are relationships between 



different values of q. For example, e x = 



-qx~ 



1 These properties reflect the possible extensivity of the nonadditive entropy S q in the presence 
of special correlations [T71 l2il l26l l28l [29] . 
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The following lemma establishes a general relationship (which contains the 
previous one as a particular case) between different g-Gaussians. 

Lemma 2.1. For any real q\, fix > and 5 > there exist uniquely determined 
92 = 3a (li j an d 02 — fciS, Pi), such that 

Moreover, q<i = S (5 — 1 + qi), 02 = <%■ 

Proof. Let q\ E R 1 , 0i > and S > be any fixed real numbers. For the equation 

(1 - (f - qi )l3 lX 2 )^ = (l-(l- q2)P2X 2 )^ 

to be an identity, it is needed (1 — qi)0i = (I — f?2 )/?2 , 1 — <?i = 5(1 — q2). These 
equations have a unique solution q2 = <5 _1 (<5 — f + q±), 02 = □ 

The set of all g-Gaussians with a positive constant factor will be denoted by 
G q , i.e., 

G q = {bG q ((3,x) :b>0,P>0}. 

2.4. g-Fourier transform 

From now on we assume that 1 < q < 3. For these values of q we introduce the 
q-Fourier transform F q , an operator which coincides with the Fourier transform 
if q = 1. Note that the q- Fourier transform is defined on the basis of the q- 
product and the ^-exponential, and, in contrast to the usual Fourier transform, is 
a nonlinear transform for q € (1,3). Let / be a non-negative measurable function 
with supp f C R. The g-Fourier transform for / is defined by the formula 

F q [f](0= I e l q ^® q f(x)dx, (2.3) 

J supp f 

where the integral is understood in the Lebesgue sense. For discrete functions 
fk, k 6 Z = {0, ±1, ...}, F q is defined as 

F q [f}(0= E e^®,/*, (2-4) 

where Zf — {k e Z : fk ^ 0}. 

In what follows we use the same notation in both cases. We also call (|2.3p or 
(|2.4[) the q- characteristic function of a given random variable X with an associated 
density f(x), using the notations F q (X) or F q (f) equivalently. The following lemma 
establishes the expression of the g-Fourier transform in terms of the standard 
product, instead of the g-product. 

Lemma 2.2. The q-Fourier transform can be written in the form 

/>oo 

F q [f](0= f(x)e^f^ q - 1 dx. (2.5) 
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Proof. For x G supp f we have 



? ®? fix) = [1 + (1 - q)ix£ + [fix)] 1 -" - 1}^ 



f(x)[l + 0.-q)ixC[f{x)]"- 1 ]^. (2.6) 
Integrating both sides of Eq. (|2.6[) we obtain l|2.5[) . □ 

Analogously, for a discrete k £ Z, (|2.4p can be represented as 
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fce,z 



Corollary 2.3. The q-Fourier transform exists for any nonnegative f G L\(R). 
Moreover, \F q [fm\<\\f\\ Ll -E 



Proof. This is a simple implication of Lemma 12.21 and of the asymptotics of e q x for 
large |x| mentioned above. □ 

Corollary 2.4. Assume f(x) > 0, x G R and F q [f](£) = for all £ € R. Then 
f(x) = for almost all x G R. 

Lemma 2.5. Let 1 < q < 3. For i/ie q-Fourier transform of a q-Gaussian, the 
following formula holds: 



i[G q W;xM) = (e q 4,J ° 9 ) . (2.7) 
Proof. Denote a — ^ an d write 

/oo 
(ae-^ 2 )® g (e^)d* 
-oo 

using the property e q +y = e x q ® q e v q of the (/-exponential, in the form 

e -0- 2 +ia"-^ dx = a ^ > ~^^ dx = 

-oo J — oo 



ia^- 1 ^ \2 a 2(g-l) g 2 



Gq ^q &q dx. 



The substitution y = ~ ^Vff y^ e ^ s ^ ne equation 



F q [ae-^](0 = ^l e~/® q e q ^ dy , 

VP 00 + 27/ 

> Q-l 

where 77 = q^/T' Moreover, the Cauchy theorem on integrals over closed curves is 
applicable because of at least a power-law decay of q-exponential for any 1 < q < 3. 



OO+ir/ _ a 2(<j-l) 5 2 



2 Here, and elsewhere, H/Hij = f-j, f(x)dx, L\ being the space of absolutely integrable functions. 
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By using it, we can transfer the integration from R + if] to R. Hence, applying 
again Lemma 12.21 we have 



og(£-i) ,2 / ° 2< ?,T 1) ; 2 \ 

F 9 [G 9 C9; a! )](0 = - 2 -^ y e, V 1 dy 

Simplifying the last expression, we arrive to Eq. (|2.7[) . □ 

Introduce the function z{s) = for s G (— oo, 3), and denote its inverse 
t G (— l,oo). It can be easily verified that z(-jtjy) = 7 and z(i) = z -l^ ■ 
Let (71 = z(q) and g_i = z^ 1 (q) . It follows from the mentioned properties of z(q) 
that 

z(-) = - and z(-) = — . (2.8) 

The function z(s) also possess the following two important properties 

z(s)z(2-s) = l and z(2 - s) + z _1 (s) = 2. (2.9) 
It follows from these properties that g_i + — = 2 . 

Corollary 2.6. For q-Gaussians the following q-Fourier transforms hold 

F q [G q (0;x)}(O = e-f*Mt\ qj = z(q), 1 < q < 3; (2.10) 

[G,^ (ftaOKO = e-^(«- 1 ^, g-i = 1 < 9 < 3, (2.11) 

1 

where /3*(s) = g~^~^i^T) T 07 "; ??Jore symmetrically, (3/ 2 ~ s f3^ 2 ~ s = K(s) with 



K{s) = fen 3 ; < J£T(s) < 1 /or s < 2, wii/i lim s ^_oc = #(2) = Oj. 

Remark 2.7. Note that /3»(s) > if s < 3 . 
Corollary 2.8. 27ie following mappings 

F q :Q q ~> Q qi , qi = z(q), 1 < q < 3, 

F 9 _! : -> G q , q-i = z^ 1 {q) 1 1 < q < 3, 
ZioZd and i/iey are injective. 

Corollary 2.9. There exist the following inverse q-Fourier transforms 
F' 1 : G qi -> G q , qi = z(q), 1 < q < 3, 

F q _\ : G q -> , g_! - ^(g), 1 < g < 3. 
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F M : Gx. -► Gl, qi = z(q), 1 < q < 3, 

91 <?1 1 



Lemma 2.10. The following mappings 

Fi:Gl^ G^, q-i = Z~ l {q), 1 < q < 3. 

1 1 9-1 

hold. 

Proof. The assertion of this lemma follows from Corollary 12 .81 taking into account 
the properties (|2.8p . □ 



Let us introduce the sequence q n = z n (q) = z(z n -i(q)),n = 1,2,..., with 
a given q = zo(q), q < 3. We can extend the sequence q n for negative integers 
n = —1, —2, ... as well putting g_ n = z~ n (q) — z~ 1 (zi- n (q)), n — 1,2, ... . It is not 
hard to verify that 

2q + n(l-q) 2(q - 1) 

^« = — n — r = 1 + o — i — tt> ^ = o,±i,±2,..., (2.12 

2 + n(l — 9) 2 — n(<7 — 1) 

which, for q ^ 1, can be rewritten as = yz - + Note that q n = 1 for 

all n = 0,±1,±2,..., if q = 1 and lim n ^,± oa z n (q) = 1 for all q ^ 1. It follows 
from (|2.12[) that q„ > 1 for all n < 2/(q — 1), q € (1,3). Moreover, obviously, 
2/(q— 1) > 1, if g € (1, 3), which implies gi > 1. Generalizing what has been said, 
we can conclude that the condition q n > 1 guarantees qu > 1, k = n — 1, n, n + 1, 
for three consequent members of the sequence (|2.12p . 

Let us note also that the definition of the sequence q n can be given through 
the following series of mappings. 

Definition 2.11. 

... A q_ 2 A q_x A q = q A qi A q 2 A ... (2.13) 
_-i ,-i _-i _-i _-i 

Z Z Z Z Z Z / f~ a\ 

... <- q_ 2 <— 9-1 <- 9o = 9 <- 9i <- ?2 <— ■•• (2-14) 
Furthermore, we set, for k — 1, 2, ... and n = 0, ±1, ... , 



and 



F~ k = F^ 1 0...0F- 1 



Additionally, for fc = we let = /• Summarizing the above mentioned 

relationships related to z n (q), we obtain the following assertions. 

Lemma 2.12. There holds the following duality relations 

q n - 1 + — = 2,n<EZ. (2.15) 

9n+l 

Proof. Making use the properties (|2.9[) . we obtain 

1 _ 1 



q n -l = z 1 (q„) = 2-z{2- q n ) = 2- — — =1 . 

z(q n ) q n+ i 



□ 



g-CLT Consistent with Nonextensive Statistical Mechanics 



9 



Lemma 2.13. The following mappings hold: 

k, n G Z: 



1qk + n I 



hm F« n g qn =g 



where Q is the set of classical Gaussians. 

Lemma 2.14. The following mappings hold: 

F< ^, 3 r F -!t 2 r F ^. 2 r £k r ^1 r ^? to i r\ 

■■■ yq-2 — * — *• ^9 — * — * yq 2 — * •■■ 

•■• <~ Sg_ 2 <~ <- £<j ^ Qqi ^ £ 92 <^ ■■• (2.17) 

Note that right sides of sequences in (|2.16p and (|2.17p cut off for n > 2/(q— 1). 



3. ^-GENERALIZATION OF THE CENTRAL LIMIT THEOREM 
3.1. (/-independent random variables 

In this section we establish a g-generalization of the classical CLT (see, e.g. [2"2"ll2"5] ) 
for independent identically distributed random variables with a finite variance. 
First we introduce some notions necessary to formulate the corresponding results. 
Let (fi, J 7 , P) be a probability space and X be a random variable defined on it 
with a density / e L q (R), v q {f) — ||/||* = f^ x ' oo [f(x)] q dx < 00. Introduce the 
density 

fq{x) = 7J\~i 

which is commonly referred to as the escort density [30 . Further, introduce for X 
the notions of q-mean, /i q = fJ- q (X) — JL^ xf q {x) dx, and q-variance a q = cr q (X — 

f-q) — ( x ~~ Mg) 2 ' fq( x ) dx, subject to the integrals used in these definitions to 
converge. 

The formulas below can be verified directly. 

Lemma 3.1. The following formulas hold true 

1. pL q {aX) = afi q (X); 

2. fiqiX-^iX)) =0; 

3. <J 2 q {aX) = a 2 a 2 (X). 

Further, we introduce the notions of q-independence and q- convergence. 

Definition 3.2. Two random variables Xi and Y\ are said 

1. q- independent of first type if 

F q [X + Y](£) = F g [X](£) ® q F q [Y}(£), (3.1) 

2. q- independent of second type if 

F q _ 1 [X + Y](Z) = F q _,[X}(0® g F q _ 1 [Y}(0, q = z(q-i), (3.2) 
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3. ^-independent of third type if 

F q _ 1 [X + Y}(0 = F g [X](0® q F q [Y](t;) ) q=z(q. 1 ), (3.3) 
where X = X 1 - l i q (X 1 ), Y = Y 1 - N (Y 1 ). 

All three types of g-independence generalize the classic notion of indepen- 
dence. Namely, for q = 1 the conditions (|3.1[) - (|3.3p turn into the well known 
relation 

F[fx * fy] = F[f x ] ■ F[f Y ] 
between the convolution (noted *) of two densities and the multiplication of their 
(classical) Fourier images, and holds for independent X and Y, If q ^ 1, then 
g-independence of a certain type describes a specific correlation. The relations 
(|3.1|) - (|3.3p can be rewritten in terms of densities. Let fx and fy be densities of 
X and Y respectively, and let fx+Y be the density of X + Y. Then, for instance, 
the g-independence of second type takes the form 




4-! 8,^ fx+ Y (x)dx = F q _, lf x }(0 ® q F q _ x [f Y }(0- (3.4) 



Consider an example of ^-independence of second type. Let random variables X 
and Y have g-Gaussian densities G q _ 1 (f3i,x) and G q _ 1 {(32,x) respectively. Denote 
7j = — 2 - q 3 ^ q 2( q _ 1 -i) , j = 1) 2. If X + Y is distributed according to the density 

i 

/ 3 _ \ 2-0-1 

G q _ 1 (8, x), where 6 = 2( q _ 1 -i) ) , then f|3.4[> is satisfied. Hence, X 

\8(7i+72)Cg_ 1 / 

and Y are g-independent of second type. 

The reader can easily modify the definition of g-independence to the more 
general case of g^-independence. For example, relation (|3.2[) in the case of g-- 
independence takes the form 

F qh _ 1 [X + Y](0=F qk _ 1 [X](0® qk F llk _ 1 [Y](0, g, = z(g fe _ 1 ). (3.5) 

Definition 3.3. Let X\. X2, Xn, ... be a sequence of identically distributed ran- 
dom variables. Denote Yjy = X\ + ... + Xjy. By definition, Xn, N = 1, 2, ... is said 
to be g^-independent (or gfc-i.i.d.) of first type if, for all N = 2, 3, the relations 

F qk [Y N - AT/ifcKO = F^Xt - n k ](£) ® qk ... ® qk F qk [X N - fi k }(0, Mfc = NAXi), 

(3.6) 

hold. 

Analogously, gfc-independence of second and third types can be defined for 
sequences of identically distributed random variables. 

Remark 3.4. For k = — 1 it follows from this definition the g-independence of a 
sequence of random variables, namely 

Fq-i [Yn - #M-i](0 = [*i ~ ®? - ®« F q -! [Xn ~ M-i](0, W = 2, 3, ... 

(3.7) 
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Example. Assume that Xn, N = 1, 2, is the sequence of identically distributed 

random variables with the associated density G q _ 1 (/?, x). Further, assume the sums 

X\ + ... + X^, N = 2,3, are distributed according to the density G q _ 1 (a, x), 
i 

where a = N 2 ~ q - 1 (3. Then the sequence Xn satifies the relation (|3.7|) for all 
N = 2, 3, thus being a g-independent identically distributed sequence of random 
variables. 

Definition 3.5. A sequence of random variables Xn, N = 1,2,..., is said to be 
g-convergent to a random variable X^ if limjv^oo F q [XN\{C) — Fq[Xoo)(0 locally 
uniformly in £. 

Evidently, this definition is equivalent to the weak convergence of random 
variables if q = 1. For 5 / 1 denote by W q the set of continuous functions <f> 
satisfying the condition \(f)(x)\ < C(l + |a:|) _ 9 =T , x € R. 

Definition 3.6. A sequence of random variables Xn with the density Jn is called 
weakly q-convergent to a random variable X^ with the density / if J R fN(x)dm q — > 
J R d f(x)dm q for arbitrary measure m q defined as dm q (x) = (f> q (x)dx, where <f> q £ 



where D N {q) = (\/ Nh' 2q -iO'2q-i)~ , N = 1,2, in the sense of Definition 13.51 
when N — > 00. Namely, the question we are interested in: Is there a q-normal 
distribution that attracts the sequence Zn ? For q = 1 the answer is well known 
and it is the content of the classical central limit theorem. 

3.2. Main results 

The formulation of a generalization of the central limit theorem consistent with 
nonextensive statistical mechanics depends on the type of g-independence. We 
prove the q-generalization of the central limit theorem under the condition of first 
type of g-independence. 

Theorem 1. Assume a sequence qk, k € Z, is given as i2.12\) with q^ € [1, 2). 
Let Xi, ...,Xn, ■■■ be a sequence of qk -independent (for a fixed k) of first type and 
identically distributed random variables with a finite q^-mean fi qk and a finite 
second {2q^ — l)-moment cr| 9fc _ 1 . 

Then Zn = Xl+ B^{ qk ) N ^ qk is qu-convergent to a q^-i-normal distribution 
as N — > 00. Moreover, the corresponding attractor is G qk _ 1 ([3k', x) with 




'N 




(3.8) 
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The proof of this theorem follows from Theorem 2 proved below and Lemma [2. 141 
Theorem 2 represents one element (fc = 0) in the series of assertions contained in 
Theorem 1. 

Theorem 2. Assume 1 < q < 2. Let X±, X^ 1 ... be a sequence of q- 
independent of first type and identically distributed random variables with a finite 
q-mean fi q and a finite second (2q — l)-moment o~ 2 q _ 1 . 

Then Zff = Xl+ ' -j^—^^Jh { s q-convergent to a q_i~normal distribution as 
N — > co. The corresponding q~\-Gaussian is G q _ 1 (f3; x), with 

f 3-g-i ykzi 



Proof. Let / be the density associated with X\ — \i q . First we evaluate F q {X\ — 
fiq) = F q (f(x)). Using Lemma [2~2l we have 

/oo />oo 
^® q f{x)dx= / f(x)e q x ^^ 1 dx. (3.9) 
-oo J — oo 

Making use of the asymptotic expansion e q = 1 + x + ^x 2 + o(x 2 ), x — > 0, we can 
rewrite the right hand side of (|3.9|) in the form 

nm = 

f{x) (l + ix^fix)}"- 1 - ^eifix)] 2 ^ + 0(x 2 e[f(x)] 2 ^)) dx = 

1 + - \i 2 °\ q -i^ q -i + o(C 2 ), £ -> 0. (3.10) 

In accordance with the condition of the theorem and the relation (2) in Lemma 13. 11 
fi q = . Denote Y, = D N (q)~ 1 (X j - = 1,2, .... Then Z w = Yj + ... + Y N . 

Further, it is readily seen that, for a given random variable X and real a > 0, 
there holds F q [aX](£) = F q [X]{a 2 - q C). It follows from this relation that F q (Yi) = 

F„\f \ ( . I ) . Moreover, it follows from the (/-independence of X\, X%, ... 

and the associativity of the q-product that 

F q [Z N ](0 =F q [f](-=J= )®,...®,F 9 [/](-_=L ) (TV factors). 

Nv 2 q-lO-2q-l V Nu 2q-1^2q-l 

(3.11) 

Hence, making use of properties of the q-logarithm, from (|3.1ip we obtain 
\n q F q [Z N }(0 = Nln q F q [f}( * ) = Nln q (l - § £ + (^)) = 



|C 2 + o(l), iV^oo, (3.12) 



locally uniformly by £. 
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Consequently, locally uniformly by £, 

lim F q (Z N ) = e q if . (3.13) 

N— »oo 

Thus, Zn is g-convergent to the random variable Z whose q-Fourier transform is 

e q ^ eG q . 



In accordance with Corollary 12.61 for q and some (3 there exists a density 
G q _ 1 (f3; x), q — z(g_i), such that F q _ 1 (G q _ 1 (/?; x)) — e q '" q ^ 2 ^ . Let us now find (3. 
It follows from Corollarv l2.6l (see (|2.11|l ) that /3* (q~i) = q/2. Solving this equation 
with respect to (3 we obtain 



y 9- 

The explicit form of the corresponding g_i-Gaussian reads as 



2C s a - 1 ^Ms) 



■2-,-: 



G q _ ± (P; x) = C- 1 ( _ fl v _ 1 m | e s J ,s = q^. (3.15) 



□ 



Analogously, the q-CLT can be proved for g^-i.i.d. of the second and third 
types. The formulations of the corresponding theorems are given below. The reader 
can readily verify their validity through comparison with the proof of Theorem 1. 

Theorem 3. Assume a sequence qk,k S Z, is given as \2.12\) with q^ £ [1,2). 
Let X±, Xjy, ... be a sequence of qk-independent (for a fixed k) of second type 
and identically distributed random variables with a finite qu-i-mean (J-q h _ 1 and a 
finite second (2qk-i — l)-moment . 

Then Zjsi = 1+ p N ( qk t / '""^ * s 3k— l - convergent to a qk~\-normal distri- 
bution as N — > co. The parameter (3k of the corresponding attractor G qk _ 1 {(3k]x) 
is 

p k =( 3 ~%:\ 2 )^ bT . (3.i6) 

Uq^C^ 1 2 ^ 

Theorem 4. Assume a sequence qk, k £ Z, is given as V2.12]) with qt £ [1, 2). 
Let Xi, ...,-Xjv, ... be a sequence of qk-independent (for a fixed k) of third type 
and identically distributed random variables with a finite qk-mean fi qk and a finite 
second (2qk — l)-moment cr| gfc _ 1 . 

Then Zn = Xl+ -Q^( qk ) N ^ qk is qk-i~ convergent to a qk-i-normal distribu- 
tion as N — ► oo. Moreover, the corresponding attractor G qk _ 1 ([3k', x) in this case 
is the same as in the Theorem 1 with (3k given in H3.8\) . 

Obviously, = 1 if and only if q = 1. This fact yields the following corol- 
lary. 
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Corollary 3.7. Let X±, Xn, ... be a given sequence of qk -independent (of any 
type) random variables satisfying the corresponding conditions of Theorems 1-3. 
Then the attractor of Zn is a q^-normal distribution if and only if qu — 1, that is, 
in the classic case. 

4. CONCLUSION 

In the present paper we studied (/-generalizations of the classic central limit the- 
orem adapted to nonextensive statistical mechanics, depending on three types of 
correlation. Interrelation showing the dependence of a type of convergence to a 
q-Gaussian on a type of correlation is summarized in Table 1. 

Table 1. Interrelation between the type of correlation, condi- 
tions for q-mean and q- variance, type of convergence, and param- 
eter of the attractor 



Corr. type 


Conditions 


Convergence 


Gaussian parameter 


1st type 


Hq h < 00, CT2q fc _l < oo 


q k — conv. 


a I 3-<j fe _i \ 2-q k _ 1 


2nd type 


IHk-l < °°> 
<729*_1-1 < 00 


qk-i — conv. 


ft = (- ^Ix-O '^ k - 1 


3rd type 


jl qk < OO, (?2q k -l < °0 


- conv. 


a ( 3-g fc _ x \ 2-9fc-i 

Wc,^ 1 1 



In all three cases the corresponding attractor is distributed according to a 
<7fc_i-normal distribution, where qk—i = £ (<7fe)- We have noticed that qk ^ qk— l 
if qk ^ 1. In the classic case both qk = 1 and z^ 1 (qk) = 1, so that the corre- 
sponding Gaussians do not differ. So, Corollary 13.71 notes that such duality is a 
specific feature of the nonextensive statistical theory, which comes from specific 
correlations of Xj, which cause in turn nonextensivity of the phenomenon under 
study. 

Now let us briefly discuss the scaling rate, important notion in diffusion 
theory. We recall that the standard Gaussian evolved in time can be calculated 

G(t,x) = F-He - ***] = -^=e~i, t > 0. 

It follows immediately that the mean squared displacement is related to time 
with the scaling rate x 2 ~ t s with 5=1. Applying the same technique in the 
case of the g-theory, in all three types of correlations, using (|3 . 1 5[) . Lemma [2.121 
and expression for (3k (see fourth column of Table 1), we obtain 6 = 1/(2 — 
qk-i) = qk+i- Hence, three consequent members of qk — 2 2+k{i~q) > 1 ^ I 1 ' 2), 
namely the triplet (qk-i,qk,Qk+i)> play an important role in the description of 
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a nonextensive phenomenon. Namely, if a correlation is given by q k , then the 
corresponding attractor is a g^-i-Gaussian and, in turn, the scaling rate is equal 
to q k +i- 

Finally, we would like to address connections of the theorems that we have 
proved with two known results. We have seen that the classical central limit theo- 
rem may in principle be generalized in various manners, each of them referring to 
correlations of specific kinds. In [3], an example of correlated model was discussed 
numerically. The correlations were introduced, in a scale-invariant manner, through 
a g-product in the space of the joint probabilities of N binary variables, with 
< q < 1@. It was numerically shown that the attractors are very close (although 
not exactly [27]) to (double-branched) Q-Gaussians, with Q = 2 — | € (— oo, 1], 
and that the model is superdiffusive [28] with 1 < S < 2). The relation Q — 2 — - 
corresponds to the particular case k = — 1 of the present Theorem 1. It comes from 

Lemma [2.121 with o_ 2 = 2 — - = 2 — which holds when k=-l. It should be 
1 ' H * qo q' 

noted the following connection between these two models, related to the behavior 
of the scaling rate S. In the model introduced in [4J, supcrdiffusion occurs with S 
monotonically decreasing from 2 to 1 when q increases from to 1 [35] . In our 
k = —1 model we have 5 = q € [1, 3), that is, monotonically increasing from 1 to 
3, when g_i increases from 1 to 2. The scaling rates in these two models behave 
as a sort of continuation of each other, through q = 1, on the interval q G [0, 3), 
and in both models we observe only superdiffusion. 

Another example is suggested by the exact stable solutions of a nonlinear 
Fokker-Planck equation in [5] . The correlation^ are introduced through aq — 2 — Q 
exponent in the spatial member of the equation (the second derivative term). The 
solutions are Q-Gaussians with Q 6 (— oo, 3), and S = 2/(3 — Q) £ [0, oo], hence 
both superdiffusion and subdiffusion can exist in addition to normal diffusion. 
This model is particularly interesting because the scaling S = 2/(3 — Q) was 
conjectured in [31j . and it was verified in various experimental and computational 
studies [32] [331 [32 [35]. 

In the particular case of Theorem 1, k — 1, we have S = q-2 = 1/(2 — q). 
This result coincides with that of the nonlinear Fokker-Planck equation mentioned 
above. Indeed, in our theorem (k = 1) we require the finitness of (2q— l)-variance. 
Denoting 2q — 1 = Q, we get S = 1/(2 — q) = 2/(3 — Q). Notice, however, that this 
example differs from the nonlinear Fokker-Planck above. Indeed, although we do 
obtain, from the finiteness of the second moment, the same expression for 5, the 
attractor is not a Q-Gaussian, but rather a ^±1 -Gaussian. 

Summarizing, the present Theorems 1-4 suggest a quite general and rich 
structure at the basis of nonextensive statistical mechanics. Moreover, they recover, 
as particular instances, central relations emerging in the above two examples. The 



3 Even though we are analyzing here the case q > 1, it is useful to compare our model with the 
one in [4]. 

4 The correlation defined in [6] is different from the g-independence of types 1-3 introduced in 
this paper. 
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structure we have presently shown might pave a deep understanding of the so-called 
g-triplet {q sen ,q re i,q s tat), where sen, re I and stat respectively stand for sensitivity 
to the initial conditions, relaxation, and stationary state [36( 1371 [3"B] in nonextensive 
statistics. This remains however as a challenge at the present stage. Another open 
question - very relevant in what concerns physical applications - refers to whether 
the ^-independence addressed here reflects a sort of asymptotic scale-invariance as 
N increases. 
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